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' The properties of neutron star matter above nuclear density are not precisely known. 

,* ^ 1 Gravitational waves emitted from binary neutron stars during their late stages of in- 

^q' spiral and merger contain imprints of the neutron-star equation of state. Measuring 

departures from the point-particle limit of the late inspiral waveform allows one to mea- 
sure properties of the equation of state via gravitational wave observations. This and a 
companion talk by J. S. Read reports a comparison of numerical waveforms from simu- 
^ ' lations of inspiraling neutron-star binaries, computed for equations of state with varying 

£Nj stiffness. We calculate the signal strength of the difference between waveforms for var- 

f\J , ious commissioned and proposed interferometric gravitational wave detectors and show 

00 ' that observations at frequencies around 1 kHz will be able to measure a compactness 

parameter and constrain the possible neutron-star equations of state. 

00 

C^s ' In simulations of the late inspiral and merger of binary neutron-star systems, one 

typically specifies an equation of state (EOS) for the matter, performs a numerical 
evolution and extracts the gravitational waveforms produced in the inspiral. In 
this talk we report work on the inverse problem: if gravitational waves from an 
inspiraling neutron-star binary are observed, can they be used to infer the bulk 



X 



properties of neutron star matter and, if so, with what accuracy? To answer this 



C3 



question, we performed a number of simulations, 1-3 using the evolution and initial 
data codes of Shibata and Uryu, while systematically varying the stiffness of a 
parameterized EOS. This parameterized EOS was previously developed in Refs. 
3,4 and is of piecewise polytropic form, p(p) = Kip Ti in a set of three intervals 
pi-i ^ P ^ Pi in rest-mass density, with the constants Ki determined by requiring 
continuity on each dividing pi and the energy density determined by the first law 
of thermodynamics. As described in Refs. 3,4, the nonpolytropic EOS of the crust 
(0 ^ p ^ po) as well as the dividing densities p\-,pi are fixed while the parameters 
{pi = p(pi), T%, T%, T3} are generally varied. In this first set of simulations we set 
Ii = = I^j = 3 and change the EOS stiffness by varying only p\, while keeping 
the Schwarzschild mass of each neutron star fixed at 1.35 Mq. The choice of EOS 
parameter varied in this work is motivated by the fact that neutron-star radius is 
closely tied to the pressure at density not far above nuclear equilibrium density. 5 
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Variation of the adiabatic exponents is the scope of a next set of simulations. 

We compared the gravitational waveforms from the simulations to point-particle 
waveforms (see, for example, Ref. 6) and calculated the signal strength of the dif- 
ference in waveforms using the sensitivity curves of commissioned and proposed 
gravitational wave detectors. We find that, as the stars approach their final plunge 
and merger, the gravitational phase accumulates more rapidly for larger values of 
pi or smaller values of the neutron-star compactness (the ratio of the neutron-star 
mass to its radius). The waveform analysis indicates that realistic EOS will result 
in waveforms that are distinguishable from point-particle inspiral at an effective 
distance (the distance to an optimally oriented and located system that would pro- 
duce an equivalent waveform amplitude) of Do = 100 Mpc or less with gravitational 
wave detectors with the sensitivity of broadband Advanced LIGO. We further esti- 
mate that observations of this sensitivity will be able to constrain p\ for a source 
at effective distance D with an accuracy of Spi/pi <~ 0.2D/D . Related estimates 
of radius measurability show that such observations can determine the radius to an 
accuracy of SR <~ 1 km D/D n . These first estimates neglect other details of internal 
structure which are expected to give smaller tidal effect corrections. This is the 
subject of work underway, which involves improving the accuracy of the estimates 
with variation of the adiabatic exponents, determination of surfaces in the equa- 
tion of state (EOS) parameter space associated with a given departure from the 
waveform of point-particle inspiral and numerical simulation of more orbits in the 
late inspiral. Also, the results mentioned above do not take into account multiple 
detectors, parameter correlation, or multiple observations. The latter possibility is 
briefly discussed below. 

In the calculations mentioned above we estimated the error fj in measuring 
an EOS parameter (such as p\ or, more precisely, a related parameter that labels 
surfaces of constant departure from point-particle inspiral) from observation of one 
event at a reference effective distance D = 100 Mpc. Here we wish to estimate the 
effect of multiple observations on the measurement accuracy. If TVj identical events, 
each with measurement uncertainty tr^, occurred at the same effective distance Dj, 
then the overall uncertainty (standard error) of the combined measurement would 
be Gil \fN~i. However, events do not occur at the same effective distance. Instead, we 
shall assume that events are homogeneously distributed in a sphere of effective ra- 
dius -D m ax — 300 Mpc. (A uniform probability distribution of events in space is also 
uniform in effective space.) We divide this sphere into / shells of effective distance 
Di = D max (i — ^)/I with i = 1, and assume that detections will only be counted 
for sources with effective distance smaller than D max . Because uncertainty scales 
linearly with effective distance, we have Gi = aoDi/D . Combining measurements 
at different distances will then result in an overall uncertainty a given by 



where Ni is the number of events in the z-th shell. Note that the JVj are random 
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variables, so a is also a random variable with some probability distribution. The 
total number of events in a given year, 

/ 

N = J2 N *> ( 2 ) 

i=l 

is itself a random variable and is Poisson-distributed around the rate of events 
1Z = (N) (average number of events per year). The probability distribution function 
V(u\TZ) of a given 1Z, and its moments (a), (a 2 ) etc., are more easily computed via 
Monte-Carlo simulation rather than analytically. However, an analytical estimate 
for the moment (er~ 2 ) can be obtained by assuming a constant density n of events 
per year per unit effective volume and converting the sums ([2]) and (jTJ) to integrals: 

(N) = Ann D 2 dD = (3) 

1 \ Dl f D — 1 . D 2 
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f -^AnnD 2 dD = ^4nnD max (4) 

Jo D 2 tr 5 

Eliminating n from eqs. © and Q yields the simple formula 

(a- 2 )-V 2 =* ^(m)-V 2 (5) 
^>o 

Our Monte-Carlo simulations confirm that, while (a), (a 2 ) 1 / 2 and other moments do 
not in general scale proportionately to 7vL~ 1//2 for fixed -D max , the moment (cj~ 2 )~ 1/ ' 2 
scales exactly as dictated by eq. ([5]). This formula indicates that, for example, three 
events randomly distributed in a sphere of effective radius -D max = 3-Do give the 
same "average" uncertainty (er -2 ) -1 / 2 = cto as one event at effective distance Dq. 
Although the above calculation was done for measurement of a single parameter, 
it can be straightforwardly generalized for more parameters, by replacing the un- 
certainty a 2 with a Fisher information matrix. As noted above, we have restricted 
consideration to neutron stars with a single fixed mass. Independent variation of 
the mass of each companion is the subject of future work. 
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